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We review our progress on the lattice calculation of low moments of both the unpolarised and polarised nucleon 
structure functions. 



1. INTRODUCTION 

Much of our knowledge about QCD and the 
structure of nucleons has been derived from 
charged current Deep Inelastic Scattering (DIS) 
experiments in which a lepton (usually an elec- 
tron or muon) radiates a virtual photon at large 
spacelike momentum —q^ = Q^. This interacts 
with a nucleon (usually a proton) and measur- 
ing the total cross section yields information on 4 
structure functions for the nucleon - -Fi , -F2 when 
polarisations are summed over and 51, (72 when 
the incident lepton beam and the proton have def- 
inite polarisations. The structure functions are 
functions of the Bjorken variable x {0 < x < 1) 
and Q^. While the unpolarised case has been 
studied experimentally for many years (since the 
pioneering discoveries at SLAG), only recently 
have experiments been reported with polarised 
beams and targets ([|l],D)- 

Theoretically, deviations from the simplest par- 
ton model (in which the hadron may be consid- 
ered as a collection of non-interacting quarks) is 
taken as strong evidence for the evidence of QCD, 
where we have interacting quarks and gluons. 
A direct theoretical calculation of the structure 
functions seems not to be possible; however using 
the Wilson Operator Product Expansion (OPE) 
we may relate moments of the structure functions 
to matrix elements of certain operators in a twist 
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or Taylor expansion in 1/Q^. For the unpolarised 
structure functions we find^, 
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where n is even starting at 2, f — u, d, s, g and 
{p, s) = 2vlP[p^^ • • -p^" - Tr], (2) 
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For the polarised structure functions, 
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dxx''gi{x, Q"^) 



iWe use {p,s|p',s'> = {2TTf2EpS(p - ■p')53^s, with = 
—m?. The Wilson coefficients E = E{ii'^ /Q^ , g{ii)) are 
known perturbatively, for example for the charged current 
we have si"' = e(9)2(l + 0(g^)), E^'^ = e^i)'^0{g'^). 
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For 51 we start with n = 0, while for 172 we begin 
with n = 2 {O^g, n = is a special case). 

While the Wilson coefficients, E, can be calcu- 
lated perturbatively, the matrix elements cannot: 
a non-perturbative method must be employed. 
In principle lattice gauge theories provide such 
a method to compute these matrix elements from 
first principles. In this talk we shall describe our 
progress towards this goal, |^,^, and try to dis- 
cuss the problems that must be overcome. We 
will not dwell on lattice details and only give re- 
sults in the form of Edinburgh (or APE) graphs 
which are plotted in terms of physical quantities. 

The moments, eq. (^, have a parton model in- 
terpretation, being powers of the fraction of the 
nucleon momentum carried by the parton 
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where q{x, /i) — q^{x, /i) -I- qi{x, fi) (and similarly 
for q{x,iJ,)) with q-^lx^fi) being the quark distri- 
butions at some scale /i and spin |. For the gluon 
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Often we re-write the distributions in terms of 
valence, v and sea, S contributions {q — u,d), 
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assuming SUf{3) symmetry. While in e(/i)iV DIS 
we can use only the charge conjugation positive 
moments (ie for even n) , in vN experiments we do 
not have this restriction and can use all moments. 
The lowest moment is particularly interesting, as 
due to conservation of the energy momentum ten- 
sor we have the sum rule 'J2q{^)^'^^ + (x)'^) = 1. 

For the polarised structure functions a similar 
interpretation holds for gi. 
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= 2 / da;x"[Ag(a;,Ai) + (-l)"Aq-(a;,^)], 
Jo 

where Aq{x,fj,) — q'j{x,fi) — qi{x,^) (and simi- 
larly for Aq{x, fi)). The suitably modified eq. (^) 
is also taken to hold. Again the lowest moment is 
of particular interest because it can be related to 
the fraction of the spin carried by the quarks in 
the nucleon. Conventionally we define A.q{fi) by 



Aq{fi) = / dx[Aq{x,fj,) + Aq{x,iJ,)], 
Jo 



and then 
2Ag = 4'''. 



(11) 
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52 contains not only a„ (the so-called Wandzura- 
Wilczek contribution to 52) but also dn-a. twist-3 
contribution. It has been argued that the n = 
moment of (72 should be 0, but it is not possible to 
check this sum rule on a lattice. Finally we note 
that (72 does not have a partonic interpretation. 

2. THE APE PLOT 

We shall now motivate the idea of an Edin- 
burgh (or APE) plot. In a lattice calculation, 
we must first Euclideanise {t — > —it) and then 
discretise the QCD action and operators. This 
introduces a lattice spacing a as an ultraviolet 
cutoff in addition to the (bare) coupling constant 
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P = G/g^ and for the Wilson formulation of lat- 
tice fermions, 1/k, which is related to the quark 
mass. (We shall only consider m„ = here.) 

We first consider measuring the light hadron 
spectrum. Numerically we measure dimensionless 
numbers representing the different masses. In a 
region where scaling is taking place, then 
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quark mass=m„/mjv,™/f 



const. 



(13) 

Thus if we plot, for example miq /nip as a func- 
tion of TO^/rrip (here we will take mN/mp = 
/{{rriTr/mp)'^)), we must be in a region where we 
see universal scaling, ie for different g values the 
results lie on the same curve. The lattice spacing 
can be found from a = mh/rmFfl^^" . 

Measuring the various masses yields the plot 
shown in Fig. Q The plot allows the whole range 
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Figure 1. The APE plot mjv/rrip against (m^/mp)^ for 
our results (filled circles, /3 = 6.0, 16^ X 32 lattice), [4]. For 
comparison we also show the results of [6], (triangles, (3 = 
6.0, 18^ X 32 lattice), [7] (inverted triangles, /3 = 6.0, 24^ x 
32 lattice), and [8] (diamonds, /3 = 5.7, 16^ X 20 lattice). 
The filled squares are the experimental result and heavy 
quark limit (1,3/2). 

of quark mass to be shown, from the heavy quark 



mass limit to the physical point. In fact from the 
lattice many quark mass worlds can be computed. 
(Unfortunately everything except the real world: 
this is because the u and d quark masses are very 
close to the chiral limit, which numerically is dif- 
ficult to realise.) Our lightest quark mass lies at 
about the mass of the strange quark (in the mid- 
dle of the plot, (m^/mp)^ « 0.5). Do we have uni- 
versal scaling? Roughly the points seem to lie in a 
band - a good sign, although the smaller coupling 
values seem to lie slightly above the larger cou- 
pling value. Are there finite size effects? Look- 
ing at the lightest quark masses we see a slight 
dropping of the value when increasing the spatial 
lattice size. This might be interpreted as a sign of 
some finite volume effects. We note also the diffi- 
culty of linear extrapolations - using our results 
alone would give a rather different result to the 
physical result. Using lighter quark mass results 
as well gives a more satisfactory extrapolation, 
1^. The curve from the heavy quark mass world 
to the physical world is certainly rather unpleas- 
ant - the curve first goes up before sinking down. 
We have also used the quenched approximation in 
which the fermion determinant is ignored (other- 
wise the computations are too costly). This in- 
troduces an uncontrolled approximation ~ for ex- 
ample it is thus not clear that the ratio in eq. ( [T^ ) 
need be the physical value. Finally as we have a 
choice of which hadron mass to use to find a, if 
we are not in the neighbourhood of the physical 
point, we will get ambiguities in the value. From 
BFl we find a « 2.0-2.5GcV"^ (using mAT, nip). 

3. LOW MOMENTS 

We shall now present our results for the low 
moments of the structure functions. The com- 
putational method is to evaluate the nucleon two 
point functions with a bilinear quark or gluon op- 
erator insertion, leading to a numerical estimate 
of the matrix elements in eqs. (H,^. There are two 
basic types: The first is quark insertion in one of 
the nucleon quark lines (quark line connected). 
As this involves a valence quark, we shall take 
this as probing the valence distributions in eq. (^ . 
In the second type the operator interacts via the 



2mjv = 0.440(4), nip = 0.321(4) at Kc = 0.15699(5). 



4 



exchange of gluons with the nucleon (for quark 
hne disconnected and gluon operators). These 
are taken as probes of the sea and gluon distri- 
butions. (Of course as we work in the quenched 
approximation there may be some distortions in 
the various distributions.) Due to technical dif- 
ficulties, such as the operator growing too large 
in comparison with the lattice spatial size and 
the mixing of the lattice operators with lower di- 
mensional operators, it has only been possible to 
measure the lowest three quark moments for the 
unpolarised valence structure functions, to give 
an estimate for the lowest gluon moment and for 
the polarised valence structure functions to com- 
pute hq, a2 and d2. (Some sea quark distributions 
also have been attempted in (|,p.) After compu- 
tation the lattice matrix elements must be renor- 
malised. For this we use at present the one-loop 
perturbation results, M. (However improvements 
are being developed, (lO|.) Finally we shall quote 
all our results at the reference scale, 
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(to avoid having to re-sum large logarithms in Z 
or the Wilson coefficients). 

Unpolarised 

In Fig. H we show the APE plot for the 3 low- 
est valence moments. Also shown are the physical 
values (taken from |p4| ) and the heavy quark mass 
limit| qi^vix) Cq^S{x-l/3), S{x) ^ 0. We see 
that the results roughly seem to lie on a smooth 
extrapolation from the heavy quark mass limit 
to the physical result. (Although the lightest 
quark mass seems to be levelling off, this should 
be checked on a larger lattice to see that there are 
no finite volume effects present.) It is to be hoped 
that smaller quark masses will continue to drift 
downwards. Certainly little sign of overshooting 
is seen (as happened in the APE mass plot). In 
Fig. ^ we show the ratios of the Uy to dy moments. 
There seem indications of a very smooth (linear?) 
behaviour with the quark mass. In Fig. ^ we show 
the lowest gluon moment. The experimental re- 
sult comes from ||lj]. (In the heavy quark mass 
limit we would expect all the nucleon momentum 
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Figure 2. The APE plot for valence results: 
with n = 2,3,4, top to bottom, in turn first for the u 
quarks (squares) and then for the d quarks (circles). (For 
{x), we have evaluated the matrix elements in two distinct 
ways, [4], methods a,b.) The experimental numbers and 
heavy quark limit are the filled symbols. 

to reside in the quarks and none in the gluons.) 
At present the signal is not so good. This is nu- 
merically a difficult calculation however and re- 
quires a large statistic. No statement at present 
can be made on the momentum sum rule. 

Polarised 

We now turn to a consideration of the polarised 
structure function moments. We first consider the 
non-singlet results for oq. The octet and triplet 
currents are given by 

3Fa - Da = Au + Ad - 2As « Auy + Ady 
Fa + Da ^ gA ^ Au - Ad ^ Auy - Ady (15) 

and thus we need only look at connected dia- 
grams. Space prevents us from showing results for 
gA, 3i^_A — Da or Fa/ Da, we refer the interested 
reader to the nice review psf . (Comparing the re- 
sults there with the experimental and the heavy 
quark mass limit we again see for gA an overshoot 
effect: we seem to have to climb up to reach the 
experimental value.) For the singlet results we 
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Figure 3. The APE plot for the ratio of the valence mo- 
ments: (x"-i)(")/(x"'i)(<*'. Method a {n = 2) are tri- 
angles. Method b (inverted triangles) and n = 3 (squares) 
are slightly displaced for clarity. 

have Au + Ad + As = AS and so we expect the 
disconnected terms to be more important. This is 
a very difficult problem, HJ^. Recent numerical 
estimates however show an encouraging trend to 
the experimental result, [ p3[ . 

Of the higher a moments we have been able to 
look at a2 . Changing to p = proton, n = neutron 
(rather than u, d), we plot 02 in Fig. |. (The 
numbers for a2 and ^2 are conveniently collated in 
|l4[.) Again we see a rather smooth behaviour be- 
tween the heavy quark mass limit and the experi- 
mental result. Indeed, perhaps due to the present 
rather large errors there, we could claim that we 
have reasonable agreement with experiment. 

Finally we turn to an estimation of d2 - the 
first non-leading twist operator in the OPE ex- 
pansion. In Fig. ^ we plot our results. While 
one might say that for n we are compatible with 
the experimental result, for p we have crass dis- 
agreement. At present we have no explanation 
for this result. (Indeed as we tend towards the 
chiral limit, the discrepancy grows.) In Fig. ^we 
compare our results with other theoretical esti- 
mates. At present the lattice result comes out a 



Figure 4. The APE plot for (x) 
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poor last. (But one should note that the most 
compatible result to the experiment appears to 
be the heavy quark mass limit.) 

Conclusions 

Although progress has been made, it is clear 
that much more remains to be done. In particu- 
lar for reliable extrapolations to the chiral limit, 
smaller quark mass simulations are necessary. 
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estimate as a diamond. The experimental value is a filled 
square and the heavy quark mass limit is a filled circle. 



